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Abstract 

Pulse detonation rocket engines (PDREs) offer potential performance improvements over con- 
ventional designs, but represent a challenging modeling task. A quasi- 1-D, finite-rate chemistry 
computational fluid dynamics model for PDREs is described and implemented. Four different 
PDRE geometries are evaluated in this work: a baseline detonation tube, a detonation tube with 
a straight extension, and a detonation tube with two types of converging-diverging (C-D) nozzles. 
The effect of extension length and C-D nozzle area ratio on the single-shot gasdynamics and per- 
formance of a PDRE is studied over a wide range of blowdown pressure ratios (1-1000). The 
results indicate that a C-D nozzle is generally more effective than a straight extension in improving 
PDRE performance, particularly at higher pressure ratios. Additionally, the results show that the 
blowdown process of the C-D nozzle systems could be beneficially cut off well before the pres- 
sure at the end-wall reaches the ambient value. The performance results are also compared to a 
steady-state rocket system using similar modeling assumptions. 
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Nomenclature 

A = pre-exponential constant in Arrhenius relation 

Cpk = specific heat at constant pressure for species k 

D = diameter 

e = total energy per unit mass 

Eact = activation energy in Arrhenius relation 

F = Troe correction factor for pressure-dependent reactions 

F c = Troe centering parameter for pressure-dependent reactions 

F = inviscid convective flux vector 

h jj? = heat of formation of species k 

H = source vector for area change 

kbr = backward reaction rate 

k e q = equilibrium constant 

kf r = forward reaction rate 

koo = high-pressure rate limit of pressure-dependent reactions 

ko = low-pressure rate limit of pressure-dependent reactions 

L = length 

Mk = molecular weight of species k 

p = ideal gas pressure, £!tLi ( pk/Mt)RT 

P r = reduced pressure for pressure-dependent reactions 

R = universal gas constant 

S{x) = cross sectional area as a function of x 

t = time from detonation initiation 

tbd = blowdown time, point in blowdown history when p ew = p a 
t zt = point in blowdown history when thrust first drops to zero 

T = gas temperature 

u = gas velocity in the x direction 

U = wave velocity in the x direction 

U = state vector 

V = volume 

W = source vector for finite-rate chemistry 

x = distance from end-wall 

Y k = mass fraction of species k, pk/p 

P = temperature exponent in Arrhenius relation 
e = nozzle area expansion ratio, S ex / S* 
y = ratio of specific heats 
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v = stoichiometric coefficient in chemical reaction 
p = total gas density, , p k 

v = volume ratio, (V dl + V se )/V d t, or (V dt + V cs + Vds)/V d t 
<p = blowdown pressure ratio, pi/p a 

Subscripts 

a = ambient condition 

c = corrected value 

cs = pertaining to the converging section of a nozzle 

C-J - Chapman-Jouguet state 

ds = pertaining to the diverging section of a nozzle 

dt = pertaining to the detonation tube 

ew = condition at the end-wall of the detonation tube 

ex = pertaining to the exit of a nozzle 

i = initial fill condition in the detonation tube 

j = reaction index 

k, l = species indices 

nr = number of reactions 

ns = number of species 

p = predicted value 

se - pertaining to a straight tube extension 

Superscripts 

m = time level index 
* = pertaining to the throat of a nozzle 

Introduction 

Pulse detonation engines (PDEs) have generated considerable research interest in recent years 
as a chemical propulsion system potentially offering improved performance and reduced complex- 
ity compared to conventional gas turbines and rocket engines. The detonative mode of combus- 
tion employed by these devices offers a theoretical thermodynamic advantage over the constant- 
pressure deflagrative combustion mode used in conventional engines. However, the unsteady 
blowdown process intrinsic to all pulse detonation devices has made realistic estimates of the actual 
propulsive performance of PDEs problematic. The recent review article by Kailasanath 1 highlights 
some of the progress that has been made in comparing the available experimental measurements 
with analytical and numerical models. 

The pulse detonation rocket engine (PDRE ) 2 has received a comparatively more limited re- 
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search interest than the pulse detonation engines intended for airbreathing use. While conceptually 
similar, the rocket application of this technology requires investigation of the performance char- 
acteristics over a wide range of ambient pressures. In conventional rocket engine designs, the 
converging-diverging (C-D) nozzle plays a critical role in converting the thermal energy of the 
combustion products in the combustion chamber to directed kinetic energy in the exhaust. 3 Thus, 
there has been considerable recent research interest in studying the use of nozzles for enhancing 
PDE performance. 4 No attempt is made here to provide an exhaustive review of nozzle research 
for PDEs, but the most recent results, relevant to the rocket application and the current study, 
should be listed. Cooper and Shepherd 5 have recently developed a partial-fill model for modeling 
the effect of straight extensions and diverging nozzles on detonation tubes at atmospheric condi- 
tions. This model makes use of the available experimental impulse data in the literature, and the 
partial-fill numerical simulations of Li and Kailasanath. 6 Numerical studies of the effects of noz- 
zles on a single detonation tube over a wider range of ambient pressures have been performed by 
Cambier 7 and Mohanraj and Merkle. 8 A study incorporating multiple tubes into a common nozzle 
was performed by Mohanraj et al. 9 Talley and Coy 10 have developed a 0-D, constant-/ model for 
studying the performance of an idealized constant-volume combustor with a nozzle over a wide 
range of pressure ratios. 

In a previous paper by the author, 11 parametric studies of the single-shot performance of a sin- 
gle, straight detonation tube were reported over a wide range of pressure ratios. A 1-D, unsteady 
method of characteristics code, employing a constant-/ assumption behind the detonation front, 
was developed for that study. Models of this type are computationally inexpensive, and may be 
used to illustrate the critical role of nozzles in the rocket PDE application. For example, a plot 
showing the specific impulse of various PDRE and steady-state rocket engine (SSRE) configura- 
tions as a function of blowdown pressure ratio, (p = pi/p a , is shown in Figure 1. The calculations 
assume a propellant gas mixture of stoichiometric H 2 -O 2 at p, = 1 atm, T t = 300 K. In each sim- 
ulation, execution stops when the pressure at the end-wall decays to the ambient value, p ew — p a . 
Note that the SSRE performance calculations employ the same constant-/ assumption used in the 
PDRE calculations. There are four performance curves in the figure: a SSRE with a sonic noz- 
zle, a SSRE equipped with an optimized C-D nozzle at each pressure ratio, a simple detonation 
tube (effectively a PDRE without a nozzle), and a PDRE fitted with a dynamically optimized C-D 
nozzle throughout the blowdown process. 

In examining the figure, a PDRE without a nozzle is perhaps best initially compared with a 
SSRE fitted with a sonic nozzle. Note that in this case the simple detonation tube provides supe- 
rior specific impulse, compared to a SSRE with a sonic nozzle, over the entire range of pressure 
ratios. Note, however, that this advantage is relatively minor when expanding to the lowest am- 
bient pressures. Moreover, in general a PDRE without a nozzle does not compare favorably to 
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a SSRE fitted with an optimized C-D nozzle at roughly <p > 1. This result is primarily due to 
the largely-choked outflow from a simple detonation tube, as there is relatively little conversion of 
thermal energy to directed kinetic energy. However, it is also possible to calculate the performance 
of a detonation tube fitted with a dynamically optimized C-D nozzle throughout the blowdown 
process. This analysis effectively assumes quasi-steady flow within the nozzle at each time step 
in the simulation. The calculations show that if such a nozzle could be fitted to a PDRE, then the 
specific impulse of the device would exceed that of a comparable SSRE over the entire range of 
pressure ratios. While such a nozzle is a considerable idealization, these simple constant-/ results 
clearly demonstrate the critical role that nozzles play in the performance of PDREs, particularly at 
the high pressure ratios typical for boost or in-space rocket applications. 

The purpose of this paper is to study the effect of straight extensions and C-D nozzles on PDRE 
gasdynamics and performance, in greater detail than can be obtained using the previous constant-/ 
method of characteristics approach. Details of the quasi- 1-D, finite-rate chemistry computational 
fluid dynamics (CFD) model developed by the author are provided first. Four different PDRE 
geometries are evaluated: a baseline detonation tube, a detonation tube with a straight extension, 
and a detonation tube with two types of converging-diverging nozzles. The effect of extension 
length and C-D nozzle area ratio on the single-shot gasdynamics, performance, and blowdown 
time of a PDRE is studied over a wide range of blowdown pressure ratios (1-1000). The overall 
effectiveness of straight extensions and converging diverging nozzles in improving performance is 
evaluated. The results are also compared to a steady-state rocket system using similar modeling 
assumptions. 


Governing Equations 

The quasi-one-dimensional Euler equations are the governing model used in this study. A 
quasi- 1-D model is at best an approximation of a pulse detonation or steady-state rocket engine. 
Real detonation waves exhibit clear multi-dimensional behavior. Additionally, flow separation in 
rocket engine nozzles is strongly dependent on multi-dimensional and viscous effects . However, 
as the primary goal of this study is to determine the effect of extensions and nozzles on pulse 
detonation rocket systems to first order, the quasi- 1-D Euler equations were deemed a reasonable 
model. The time-dependent form of these equations, including the effect of finite-rate chemistry, 
is written as follows: 

2u + ^E5 = h + w (d 

dt S(x) dx 
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where the state vector, U, and inviscid convective flux vector, F, are given by 
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and the source vectors accounting for the effect of area change, H, and for the effect of finite-rate 
chemistry, W, are 


0 


lUl 

0 


U>2 

• 

w = 

* 

0 


w ns 

£35 
5 dx 


0 

0 


0 


The total energy per unit mass, e, is given by 


and is used for the implicit evaluation of temperature by Newton-Rhapson iteration. 

In this work, the ideal gas thermodynamic fits of McBride et al. 12 are used to calculate the 
internal energy, specific heat, and Gibbs free energy for the ns = 9 species in the problem (N 2 , 02 , 
H 2 , OH, H 2 0, H, O, H0 2 , and H 2 0 2 ). 


Chemistry Model 

The source terms for finite-rate chemistry, Wk in Eq. 3, are calculated from the sum of relevant 
chemical reactions for each species. Thus, for a chemical mechanism of nr elementary reactions, 
each reaction is expressed as 

ns ns 

/= 1 /=i 

where the vjj are the stoichiometric coefficients of species l on the reactant side of reaction j, while 
the v'ij are the corresponding stoichiometric coefficients on the product side. The production term 
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for each species k may be written as 



Details of the chemical kinetics mechanism for H 2 -O 2 combustion used in this work are shown 
in Table 1. This nr = 18 reaction mechanism was developed by Petersen and Hanson 13 to model 
H 2 -O 2 ignition at the elevated pressures typical of practical high-speed propulsion systems. The 
rate coefficients for reactions without pressure dependence take on the conventional Arrhenius 
form: 

k frj = AfT* exp (-Eacj/RT) (7) 

Pressure-dependent reactions are modeled using the approach given by Kee: 14 


kfrj = k 00 [P r /(l+P r )]F 


where the reduced pressure, P r , is 


Pr = 


*o[M] 


'■00 


and the correction factor, F, is in the Troe 15 form: 


In F = 


, . lnF r + c 

1 -| In F, 

n - c?(lnP r -I- c) 


where F c is given by 


F c = (1 -a)exp(-r/r***) + aexp(-r/r*) + exp(-r*VF) 


( 8 ) 

(9) 


( 10 ) 


(ID 


and the constants c = —0.4 - 0.67 In F c , n = -0.75 - 1.27 In F c , and d = 0.14. 

All reactions in the mechanism are reversible. The reverse reaction rate for each reaction, kt, r , 
is calculated from the forward rate and the equilibrium constant: 

hr, = ^ 02 ) 

K eqj 

where the equilibrium constant, k eqj , is calculated using the standard Gibbs free energy difference 
for each reaction. 
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Numerical Method 

The governing quasi- 1-D Euler equations are solved in finite-volume form throughout the entire 
domain. The equations are solved in split form, in which a fluid convection subroutine solves Eq. 1 
without the chemistry source vector W, 


au 1 3FS 

3/ + S(*) dx 


(13) 


and a finite-rate chemistry integration subroutine solves Eq. 1 as a system of ordinary differential 
equations ignoring F and H, 



(14) 


The code utilizes the Strang 16 second-order time step splitting approach to couple the fluid 
convection and finite-rate chemistry subroutines. For each complete time step, the chemistry sub- 
routine is first called for one-half of a time step, followed by a full fluid convection time step, 
and then followed by another half-step of the chemistry routine. The time step splitting approach 
allows both high-quality fluid and chemical solvers to be developed and tested independently, and 
then joined together in relatively straightforward fashion. As described by Oran and Boris, 17 time 
step splitting works well when relatively small time steps are used. Since the global time step in 
the present code is strictly limited by the explicit Courant-Friedrichs-Lewy (CFL) condition, this 
is not thought to be a problem. 

The fluid solver used here is the explicit, 2nd-order accurate (in time and space), symmetric- 
TVD algorithm described by Yee. 18 The solver employs Roe’s approximate Riemann solver 
modified for nonequilibrium ideal gases. 19 It also incorporates suggestions by Larrouturou 20 to 
ensure species positivity. 

As with most reactive flow problems, the time integration of the chemistry mechanism in this 
work requires a stiff-ODE solver to ensure accuracy. The method used here employs a semi- 
implicit trapezoidal scheme as a predictor: 


[i - ^Arjj (Up +1 - U m ) = A tW m (15) 

where I is the identity matrix, and J is the Jacobian of W, J = 3W/3U. Newton iteration of the 
trapezoidal scheme is used as a corrector: 



(U« +1 - U™ +1 ) = IT - U£ +1 + ^Ar(Vr + W£ +1 ) 

w 


(16) 


The chemistry integration subroutine is designed to use several sub-steps, if necessary for accu- 
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racy, within a given global half-step. The sub-steps are chosen such that no species concentration is 
predicted to change by more than 1% within each sub-step. Naturally, this leads to more extensive 
use of sub-steps when the concentration within a computational cell is far from equilibrium. The 
Jacobian J is evaluated analytically once per sub-step at U m . Thus, the LU decomposition of the 
bracketed matrix terms on the left hand side of Eq. 15 can be re-used in Eq. 16. Two corrector 
iterations were performed per sub-step. 

PDRE Simulations 

The four different PDRE geometries studied in this work (Fig. 2) are highly idealized. Each is 
based on a constant-area (£><*, = 2.0 cm) detonation tube = 20 cm in length. One end of the 
detonation tube is closed, while the other is either open to the environment (geometry A), attached 
to a straight extension (geometry B) or a converging-diverging nozzle section (geometries C and 
D). Geometry A represents the baseline case. Four different straight extensions for geometry B 
were evaluated in this study: L se = 4.0 cm, 10cm, 20 cm, and 40 cm. These extensions correspond 
to ratios of the total assembled tube volume to the detonation tube volume of v = 1.2, 1.5, 2.0 
and 3.0. In both of the converging-diverging nozzle geometries (C and D), the half-angle in the 
converging section is 15°. The primary difference between the two geometries is that in C the 
length of the converging section is L cs — 0.4 cm, while in D this length is L cs = 1.2 cm. This 
results in a nozzle throat diameter for geometry C of D* = 1.786 cm, and for geometry D of 
D* = 1.357 cm. Thus, the maximum rate of outflow is more restricted in geometry D than in C, 
and in both cases is less than that of baseline geometry A. The half-angle in the diverging section of 
the nozzle is fixed at 15° in both geometries. The length of the diverging section, is a variable 
set at run-time, as detailed in the results and discussion to follow. The C-D nozzle area ratio is 
defined throughout this work as the ratio of the nozzle area at the exit plane to the throat area. 

The detonation tube is pre-filled with a gaseous propellant mixture with no initial velocity. 
Stoichiometric H 2 -O 2 at p, = 1 atm, T-, = 300 K is utilized for all calculations in this work. An 
idealized, massless diaphragm isolates the propellant mixture from the nozzle and ambient envi- 
ronment until ruptured by the detonation wave. The extension or nozzle section is initially filled 
with H 2 gas at the specified ambient pressure. The choice of the ambient gas has a nontrivial 
impact on the impulse of geometries B, C and D in single-shot simulations. Moreover, a number 
of options can be reasonably selected. Hydrogen is used here primarily due to its low molecular 
weight, which allowed a wide range of ambient conditions to be studied using the thermochemical 
model explained previously. The impact of using other possible choices will be addressed in a 
limited form in the results and discussion section. 

Detonation initiation is facilitated by calculating the equilibrium combustion conditions (com- 
puted at constant density and internal energy) for stoichiometric H 2 -O 2 at 10 atm and 3000 K, and 
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then specifying these conditions in the cell adjacent to the closed end-wall. 

Ghost cells are utilized to specify the boundary conditions in the problem. A reflection-type 
boundary condition is utilized at the closed end of the detonation tube to simulate a solid wall. The 
method of characteristics approach of Poinsott and Lele 21 is used to calculate the exit flow bound- 
ary condition. For sonic (choked) or supersonic exit flow, all exit flow properties are determined 
by the interior flow. For subsonic exit flow, the ambient pressure is specified and the method of 
characteristics is used to compute the remaining flow properties. If reverse flow is detected at the 
exit, a subsonic inflow boundary condition is calculated using the ambient H 2 gas as a constant 
pressure and enthalpy reservoir. 

T\vo additional special restrictions are imposed in the simulations. In order to simulate the 
effect of an idealized diaphragm, only the detonation tube portion of the domain (from the closed 
wall to the diaphragm location) is computed initially. A reflection-type boundary condition is 
specified at the diaphragm location until the pressure in the adjacent cell rises 1.0% above the initial 
fill value. This special restriction is subsequently removed, and the entire domain is computed. 
Additionally, there is a check performed when the exit flow is supersonic. Since in this case the 
exit boundary conditions are entirely calculated from the interior flow, there is no way for the exit 
flow to return to a subsonic condition. 8 Therefore, at each time step a check is made to determine if 
the pressure from a standing normal shock at the exit is less than the ambient pressure. If true, then 
the normal shock properties are specified in the last interior cell, and a subsonic outflow boundary 
condition computed at the exit. 

Grid Resolution 

Grid resolution is a critical issue in all finite-rate chemistry CFD models. Due to the extremely 
fast kinetics of the H 2 -O 2 system, numerical modeling of the reaction zone of undiluted H 2 -O 2 
Chapman-Jouguet (C-J) detonations requires an exceptionally fine grid spacing. The C-J detona- 
tion velocity of stoichiometric H 2 -O 2 at p\ - 1 atm, 7/ = 300 K is calculated by the NASA CEA 
thermochemical code 22 to be 2837.1 m/s. Using the calculated frozen shock properties at this ve- 
locity, and the chemical kinetics mechanism described previously, the length of the steady C-J 
reaction zone for this mixture can be estimated as ~ 0.02 mm. Thus, a grid spacing of roughly 
1 /xm would be required to resolve the reaction zone. This is not currently practical, nor is it nec- 
essary for the performance estimates of this study. The chief requirement here is that the grid 
resolution accurately calculate the C-J state and velocity. 

A grid resolution comparison of pressure and temperature profiles of a propagating C-J det- 
onation, t = 30 /xs after initiation, is shown in Fig. 3. The resolutions vary over one order of 
magnitude. As even the finest resolution here, Ax = 0.04 mm, is twice the estimated reaction zone 
size, none of the resolutions can accurately capture the ignition process behind the leading shock 
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wave. Rather, in all cases, the propagating C-J detonation is smeared over several cells, resulting 
in the pressure spike and subsequent relaxation toward the C-J state in the finer grids. All four res- 
olutions capture the detonation velocity with reasonable accuracy (~ 2836 m/s for Ax = 0.04 mm, 
~ 2840 m/s for the other resolutions), though the precision of this estimate is naturally lower for 
the coarser grids. Note, however, that in the two finer grids the pressure relaxes to the C-J value in 
a significantly shorter distance than in the coarser grids. Based on these observations, a uniform 
grid spacing of Ax = 0.10 mm was utilized for all simulations in this study. 

Thrust Calculations 

The time-dependent thrust is calculated at each time step by two different methods. One mea- 
sure of the thrust assumes a control volume tightly bounding the solid surfaces of the PDRE, and is 
determined by integrating the pressure difference across all surfaces. A second measure assumes 
a rectangular control volume encapsulating the PDRE. This measure of thrust is determined from 
the sum of the time rate of change of the internal momentum integral across the domain, the mo- 
mentum flux from the nozzle section, and the pressure difference across the control volume. In 
general, there is excellent agreement between the two thrust calculations, and the time-integrated 
impulse calculations agree to within 0.1%. In all simulations, the calculation proceeds until the 
pressure at the closed end-wall decreases to the ambient value (p w = p a ). Thus, these simulations 
should be thought of as single-shot results. 

Results and Discussion 

There are several critical factors which can impact the performance of a PDRE system. The 
blowdown pressure ratio is one of the most important parameters governing the performance of 
any rocket-type system. Additionally, the length of a straight extension, or the area ratio of a C- 
D nozzle, can have a large effect on PDRE performance. In the following section, the essential 
gasdynamics of the blowdown process in the four PDRE geometries are reviewed. This is followed 
by a parametric study of the effect of blowdown pressure ratio, extension length, and C-D nozzle 
area ratio on specific impulse and blowdown time. The relative effectiveness of straight extensions 
and C-D nozzles in improving performance is compared. A comparison of the performance of a 
PDRE fitted with an optimized C-D nozzle with an equivalent steady-state rocket engine is also 
provided. 

PDRE Blowdown Gasdynamics 

The blowdown process of a PDRE can be thought of as occurring in several different phases, 
typically differentiated by the location of important wave structures in the tube. Consider the 
simplest geometry, the baseline detonation tube (geometry A), using the standard initial conditions 
and p a = 0.1 atm (<p = 10). Pressure profiles, at several representative points in the blowdown 
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history, are shown for this case in the upper panel of Fig. 4. The detonation is initiated at the 
end-wall at t =0. The subsequent propagation of the detonation wave is clearly shown in the 
pressure profile labeled t = 60 /a s. The sharp pressure rise associated with the detonation front 
is followed by a Taylor expansion wave which acts to slow and expand the combustion products 
behind the detonation front. Well behind the detonation front, the products reach a steady pressure 
at zero forward velocity. Once the detonation reaches the exit of the tube and breaks the idealized 
diaphragm, outflow begins, and an expansion wave propagates back upstream. This expansion 
wave expands and accelerates the combustion products toward the exit, and is visible in the t — 
180 /as profile. Once this expansion reaches the end-wall, it reflects and begins propagating back 
downstream toward the exit again ( t — 300 /is curve). Eventually, the expansion wave reaches the 
tube exit, and further reduces the pressure at that location (t = 420 /as). Note that the pressure at 
the end-wall is steady at a “plateau” value until the primary expansion wave reaches the end-wall. 
Similarly, the pressure at the exit plane also rapidly achieves a quasi-steady value, which holds 
until the primary expansion wave again reaches the exit after reflection from the end-wall. 

The addition of an extension to the detonation tube (geometry B) primarily serves to increase 
the amount of time required for some of the waves to traverse the tube. Pressure profiles at repre- 
sentative times in the blowdown history of a detonation tube with a L se = 10 cm extension section 
are shown in the lower panel of Fig. 4. As above, ^ = 10. As in the baseline case, a primary ex- 
pansion wave is initiated when the detonation reaches the end of the detonation tube. Additionally, 
as there is no combustible gas mixture in the extension, the detonation wave changes into a trans- 
mitted shock wave in the extension. When the transmitted shock reaches the exit of the extension 
section, it initiates a secondary expansion wave which also travels upstream toward the end-wall. 
This leads to pressure profiles quite similar to the baseline case. However, the overall effect of the 
extension is to slow the rate of pressure relaxation at the end-wall. 

A converging-diverging nozzle has a significant effect on the pressure distribution of a PDRE. 
Pressure profiles are shown for two C-D nozzle example cases in the upper (geometry C with 
Lds — 3.6 cm) and lower panels (geometry D with Lds = 2.8 cm) of Fig. 5. In both cases (p is again 
10. As was the case previously, a primary expansion wave, initiated when the detonation reaches 
the end of the detonation tube, is a key factor in controlling the pressure relaxation at the end-wall. 
However, it is clear from examining the figures that the C-D nozzle acts to expand the combustion 
products to a much lower exit pressure than would otherwise be accomplished by the detonation 
tube alone. The primary difference between the two C-D nozzle cases is the significantly smaller 
throat in geometry D, which results in reduced mass outflow compared to geometry C. This flow 
restriction is significant enough that the primary expansion wave is led by a weak shock wave, 
clearly visible in the lower panel of the figure. 

The impulse history of all four of these example cases is shown in Fig. 6. The instantaneous 
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thrust can be thought of as the slope of these impulse curves. As discussed previously, the primary 
effect of the extension in geometry B is to reduce the rate of pressure relaxation at the end-wall. 
This results in a modestly higher thrust in the later portions of the blowdown history, compared 
to the baseline geometry A. In contrast, the C-D nozzle has a more significant effect on the thrust 
and impulse. Geometry C has a comparatively higher thrust and resultant impulse throughout the 
entire outflow portion of the blowdown history. Due to the more restricted mass outflow, and 
internal shock wave reflections, the instantaneous thrust in geometry D is generally lower than 
that of the other cases. However, although the overall blowdown time is significantly longer than 
geometry C, the final impulse is nearly equivalent. One additional important observation which 
can be made from Fig. 6 is that typically the integrated impulse will achieve a value quite close to, 
or even slightly exceeding, the final value well before the end of the simulation. This is particularly 
true of the two C-D nozzle geometries, and is observed over the entire blowdown pressure range 
studied in this work. This point will be addressed in greater detail in the following sections. 

PDRE Performance Results 

The performance trends just discussed can be seen in greater detail by studying the performance 
of the various PDRE geometries over a wide range of blowdown pressure ratios. The specific im- 
pulse of the baseline detonation tube (geometry A), and the detonation tube with several extension 
lengths (geometry B with L se = 4, 10, 20, and 40cm) is shown in Fig. 7. As is evident from the 
figure, the general performance trends of all five configurations are quite similar. There is a fairly 
significant improvement in performance (30-38%) over <p = 1-10, while decidedly less so (4-8%) 
over <p = 10-100. There is very little improvement in performance (~ 1%) over <)> = 100-1000. 
In general, there is a definite performance improvement provided by the extension section, and 
this benefit increases with extension length. The one exception to this trend occurs at </> = 2. In 
this particular case, there is no additional performance benefit for L se > 10 cm. This is due to 
the fact that the primary expansion wave, generated at the interface between the detonation tube 
and the extension section, is alone sufficient to reduce the pressure at the end-wall to the ambient 
value. In this particular case, this event occurs before the secondary expansion wave (generated 
when the transmitted shock reaches the exit of the extension) can reach the end-wall. Thus, there 
is no effect from the longer extensions at this pressure ratio. At higher blowdown pressure ratios 
(<p = 5-1000) this phenomena does not occur; the pressure decay (to ambient) provided by the 
primary expansion at the end-wall is long enough that the secondary expansion wave can reach the 
end-wall in time to affect the pressure relaxation rate. 

Both the blowdown pressure ratio, and the nozzle area expansion ratio, e, play critical roles 
in the performance of the C-D nozzle PDRE geometries. Performance results as a function of C- 
D nozzle area expansion ratio for geometries C and D (upper and lower panels, respectively) are 
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shown in Fig. 8. In each case, e, was adjusted by varying the length of the diverging section of 
the nozzle (L^). Results for 4> = 1-1000 are plotted. As is the case for conventional rocket-type 
systems, for a given pressure ratio, the performance gain provided by a C-D nozzle will generally 
increase with e up to a certain optimum point. After this point performance will decrease with 
increasing e. The dashed curve in each graph shows the fixed e which yields optimum specific 
impulse for each <p. It is important to recognize that this optimum point is specific to the C-D 
nozzle geometry and blowdown time assumption used. 

In general, there is little performance difference between the two C-D nozzle geometries. The 
performance of geometry D, for a fixed <j> and e, tends to be slightly lower than geometry C, 
primarily due to losses induced by the reflected shock waves present in the blowdown history of 
D. Note that a C-D nozzle is generally more effective at improving performance at higher pressure 
ratios than the extensions in geometry B. The specific impulse using an optimized C-D nozzle 
improves by 53-57% over <p = 1—10, and 20-21% over <p = 10-100. Even over 0 = 100-1000, 
the performance increases 10%. This is largely due to the ability of a C-D nozzle to convert the 
thermal energy of the combustion products into directed kinetic energy. This effect can also be 
seen by plotting specific impulse for all four geometries as a function of volume ratio, v (Fig. 9). <p 
is 1000 in the upper panel of the figure, and 10 in the lower panel. It is evident from the figure that 
a C-D nozzle is generally a more effective use of volume than a straight extension, particularly at 
relatively high <p. It should be noted, however, that this observation is based on the assumption that 
the volume ratio range under consideration is within the range where increased nozzle expansion 
provides improved performance. 

PDRE Blowdown Time Results 

As the cycle time of an operational PDRE system is of great practical interest, it is useful to 
examine the effect of the blowdown pressure ratio, straight extension length, and C-D nozzle area 
expansion ratio on the blowdown time of the PDRE geometries considered here. The blowdown 
time, tbd, in this work is defined as the time at which the pressure at the end-wall decays to the 
ambient value (p ew = p a ). The strong dependence of tbd on <p and extension length can be seen 
in Fig. 10. The time in that figure has been normalized by Ldt/Uc-j = 70.495/zs, which is the 
characteristic time for a C-J detonation to traverse the detonation tube. Note that at the highest 
blowdown pressure ratios, there is an almost directly proportional dependence of blowdown time 
to extension length. The obvious kink in the L se = 20 cm and L se = 40 cm curves at <j> — 2 was 
explained in the previous section. 

The blowdown time for the two converging-diverging nozzle PDRE geometries is shown in 
Fig. 11. The values have again been normalized by Ldt/Uc-j • As would be expected, tbd is con- 
sistently longer for geometry D than C, over the entire range of blowdown pressure ratios. Note 
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that these results assume the use of the optimum fixed C-D nozzle for each pressure ratio, though 
since the nozzle throat outflow is nearly always choked, tbd is nearly independent of e at a given <f>. 
One of the notable results of these simulations is that the instantaneous thrust of even an optimum 
fixed C-D nozzle is slightly negative for much of the latter portions of the blowdown time. This 
result is due to the fact that the C-D nozzle overexpands the residual combustion products when 
the pressure in the detonation tube has decayed sufficiently. Thus, it is useful to also plot the time 
at which the thrust first drops to zero in the blowdown history (t a ). It is evident from the figure that 
this point occurs considerably earlier than tbd, particularly at the highest <f>. As the instantaneous 
thrust, even while negative, is so small at this point, there is little difference in integrated impulse 
between these two time points (typically at most a loss of ~ 1 s of specific impulse over the in- 
terval). However, this observation suggests that a considerable savings in overall cycle time could 
be realized with an actual slight gain in performance. Note that as the end-wall is the only thrust 
surface for geometries A and B,t a = tbd for those cases. 

Effect of Ambient Gas on Results 

As these are single-shot simulations, it is logical to expect that the choice of the ambient gas 
in the extension and nozzle sections may influence the performance and blowdown time results. 
In an operational multi-cycle PDRE, the gas in these sections would likely be a mixture of the 
ambient atmosphere, the combustion products from the previous cycle, and any purge gases used 
in the cycle. Hydrogen is used as the ambient gas in all simulations in this work. However, in 
order to form a limited appreciation for the impact of this choice on the results, simulations using 
helium (He) and N 2 as the ambient gas are also conducted at the p a =0.1 atm (0 = 10) condition. 
In effect, due to their higher density. He, and particularly N 2 , provide a larger tamping mass for 
a given extension or nozzle volume than H 2 does (see Ref. 5 for a discussion of this point). In 
general, the performance results using He are only slightly higher (0-1 s of specific impulse) than 
those using H 2 . However, the use of N 2 as the ambient gas yields more significant deviations from 
the H 2 reference calculations. The greatly reduced sound speed of N 2 compared to H 2 results in 
a much stronger transmitted shock in the extension or C-D nozzle sections. In the case of a C-D 
nozzle, the resultant high pressure behind this shock causes additional impulse to be transmitted 
to the diverging section of a C-D nozzle (typically 3-4 s of additional specific impulse for this 
condition). Blowdown time was not affected. In the case of a straight extension, the stronger 
transmitted shock significantly reduces the pressure variation across the primary expansion wave. 
This causes a comparatively higher pressure at the end-wall before the arrival of the secondary 
expansion wave (from the exit). Thus, both impulse and blowdown time are increased (from ~ 
0.3% additional impulse and blowdown time at L se = 4 cm, to 12% additional impulse and 18% 
additional blowdown time at L se = 40 cm). 
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Equivalent SSRE Model and Performance Comparison 

It is instructive to compare the specific impulse of a PDRE fitted with a C-D nozzle with a 
SSRE under equivalent modeling assumptions. While frozen and equilibrium rocket performance 
calculations can be readily obtained from the CEA code, these results are not directly comparable 
to the finite-rate chemistry model used in the current PDRE code. The primary concern is the 
tendency for chemistry to slow down in real nozzle systems as the temperature and pressure are 
reduced in the expansion process. Thus, it is best to compare the finite-rate PDRE systems with 
a finite-rate SSRE model. This SSRE CFD model is heavily derived from the PDRE code. The 
Ldt — 20 cm detonation tube is replaced with a thrust chamber 0.1 cm in length, and 4.0 cm in 
diameter (Fig. 12). This thrust chamber then converges at a constant 15° angle, over a length of 
4.13 cm, to a throat D* — 1.787 cm in diameter. After the throat, the nozzle diverges at 15° to a 
length specified by the user. Thus, the throat diameter and nozzle expansion rate of the SSRE model 
closely matches that of PDRE geometry C. The equilibrium temperature, pressure and composition 
of stoichiometric H 2 -O 2 , burned at constant pressure and enthalpy from p,- = 1 atm and 7) = 300 K, 
form a reservoir for a subsonic inflow boundary condition into the domain. Thus, this analysis 
assumes that the propellant feed system of the SSRE model is the same as the PDRE. 

Similar to the PDRE nozzle optimization study, the finite-rate SSRE CFD model is run at 
various e until the optimum specific impulse is obtained for a given <j> . However, as the solution 
converges to steady-state in this model, flowfield information from the highest pressure ratio case 
can be used to guide area selection at the lower pressure ratios. In each case, the SSRE model is run 
until the solution converges. In general, the specific impulse and optimum expansion ratio using 
finite-rate chemistry are slightly larger than that for frozen chemistry. Additionally, if the chemistry 
is frozen in the SSRE CFD model, there is excellent agreement (to within 0.3% in specific impulse) 
with the frozen-chemistry predictions of CEA. 

The mixture-based specific impulse for the baseline detonation tube (PDRE geometry A), 
PDRE geometry C fitted with an optimized fixed C-D nozzle (at each blowdown pressure ra- 
tio), and the equivalent optimized SSRE system is plotted for 0 =1-1000 in Fig. 13. As would 
be expected from the results discussed previously, the relative gain from the C-D nozzle systems 
becomes more pronounced at higher pressure ratios. It is interesting to note that that both the 
baseline and C-D nozzle PDRE systems outperform a SSRE at <p < 7. This is due to the fact that 
the relative gain from the extra pressurization provided by the detonation process is greatest in this 
regime. Thus, the PDRE may have considerable theoretical potential for rocket-type applications 
when the pressure of the ambient environment is high. Additionally, at higher <p, a PDRE with 
an optimized fixed C-D nozzle still has a greater specific impulse than an equivalent SSRE. These 
results indicate that a PDRE equipped with a fixed C-D nozzle, though obviously not as efficient 
at generating thrust as a (theoretical) dynamically adaptive nozzle, can still yield single-shot per- 
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formance superior to an equivalent SSRE over a wide range of pressure ratios. Note, however, that 
this performance gain becomes relatively smaller at higher pressure ratios. 

Conclusions 

A quasi- 1-D, finite-rate chemistry CFD model for studying PDRE gasdynamics and perfor- 
mance is described and implemented. The single-shot performance and blowdown time character- 
istics of four different PDRE geometries are studied, over a range of blowdown pressure ratios from 
1-1000. The pressure at the end-wall is allowed to decay to the ambient value for all simulations 
in this work. The results show that both straight extensions and converging-diverging nozzles can 
provide valuable improvements to the performance of a baseline detonation tube. However, opti- 
mized C-D nozzles are generally more effective in increasing performance than straight extensions, 
particularly at higher pressure ratios. Studies of the blowdown time of the PDRE systems show a 
strong dependence on pressure ratio for all geometries. Straight extension length will also gener- 
ally increase blowdown time. Examination of the blowdown history of the optimized C-D nozzle 
geometries shows that the blowdown of these systems could be cut off well before the end-wall 
pressure reaches the ambient value, with an actual slight gain in performance. The performance 
of a PDRE fitted with an optimized fixed C-D nozzle is also compared with a SSRE using similar 
modeling assumptions. The results show that while a PDRE enjoys the largest relative performance 
gain over a SSRE at low pressure ratios, there is still a noticeable relative advantage for the PDRE 
at blowdown pressure ratios as high as 1000. 
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List of Table Captions 

Table 1 : Chemical kinetics mechanism for H2-O2 combustion, from Petersen and Hanson . 13 
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Number 

Reaction" 


A b 


fi b 

Eacf 

1 

O + H 2 ^ H + OH 


5.00 x 

10 4 

2.70 

6290 

2 C 

h + o 2 im^ ho 2 + M 


2.80 x 

10 18 

-0.90 

0 

3 

H + O 2 + O 2 HO 2 4- O 2 


3.00 x 

10 20 

-1.70 

0 

4 

H + 0 2 + H 2 0 ^ H0 2 +H 2 0 


9.38 x 

10 18 

-0.80 

0 

5 

H + O 2 + N 2 ^ho 2 +n 2 


2.60 x 

10 19 

- 1.20 

0 

6 

H + 0 2 ^O + OH 


8.30 x 

10 13 

0.00 

14413 

7 

H + HO 2 ^o 2 + h 2 


2.80 x 

10 13 

0.00 

1068 

8 

H + H0 2 ^ OH + OH 


1.34 x 

10 14 

0.00 

635 

9 

H + H 2 02 ^H 02 + H 2 


1.21 x 

10 7 

2.00 

5200 

10 

OH + H 2 ^ H 2 0 + H 


2.16 x 

10 8 

1.50 

5200 

ll* e 

OH + 0H + M^H 2 0 2 + M 

koo 

7.40 x 

10 13 

-0.40 

0 



ko 

2.30 x 

10 18 

-0.90 

-1700 

12 

OH + H0 2 ^ 0 2 + H 2 0 


2.90 x 

10 13 

0.00 

-500 

13^ 

OH + H 2 0 2 ^ H0 2 + H 2 0 

k a 

1.75 x 

10 12 

0.00 

320 



kb 

5.80 x 

10 14 

0.00 

9560 

14* 

H0 2 + HO 2 ^02 + H 2 O 2 

k c 

1.30 x 

10 11 

0.00 

-1630 



k d 

4.20 x 

10 14 

0.00 

12000 

15* 

o + o + m^o 2 + m 


1.20 x 

10 17 

- 1.00 

0 

I6 d 

O + H + M^OH+M 


5.00 x 

10 17 

- 1.00 

0 

17'' 

H + OH + M ^ H 2 0 + M 


2.20 x 

10 22 

- 2.00 

0 

18' 

h + h + m^h 2 + m 


1.00 x 

10 18 

- 1.00 

0 


Note: species are N 2 , O 2 , H 2 , H 2 O, OH, H, O, HO 2 , H 2 O 2 . 


“All reactions are reversible. 

b kf r = AT & exp(— E act /RT); units are in cal, mol, cm 3 , and s. 

C M does not include O 2 , H 2 O, or N 2 ; all collision efficiencies = 1.0. 
^Collision efficiencies for M; H 2 = 2.0, H 2 O = 6.0, all others = 1.0. 
*Troe parameters: a = 0.7346, T*** = 94, T* = 1756, T** = 5182. 
-^Rate coefficient is non-Arrhenius: &13 = k a + kb 
*Rate coefficient is non-Arrhenius: £14 = k c + kd 
^Collision efficiencies for M; H 2 = 2.4, H 2 O = 15.4, all others = 1.0. 
'Collision efficiencies for M; H 2 = 0.73, H 2 O = 3.65, all others = 1.0, 
' Collision efficiencies for M; H 2 = 1.7, H 2 O = 7.0, all others = 1.0. 
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List of Figure Captions 


Figure 1: Constant- y performance comparison of a single-shot PDRE with a conventional 

SSRE equipped with both sonic and optimized converging-diverging nozzles. Propellant mix- 
ture: stoichiometric H 2 -O 2 at pi = 1 atm, J,- = 300 K. 

Figure 2: Schematic of PDRE geometries studied. 

Figure 3: Grid resolution comparison showing pressure (upper panel) and temperature 

(lower panel) profiles of a propagating C- J detonation wave at t = 30 fis after initiation. 

Figure 4: Pressure profiles for PDRE geometry A (upper panel) and B with L se = 10 cm 

extension (lower panel) at representative points in the blowdown history. Propellant mixture: 
stoichiometric H 2 -O 2 at pi = 1 atm, Ti = 300 K ,p a = 0.1 atm (0 = 10). 

Figure 5: Pressure profiles for PDRE geometry C with L* = 3.6 cm (upper panel), and D 

with Lds — 2.8 cm (lower panel), at representative points in the blowdown history. Propellant 
mixture: stoichiometric H 2 -O 2 atp, = 1 atm, T, = 300 K. p a = 0.1 atm (0 = 10). 

Figure 6: Impulse history for example PDRE geometries. Propellant mixture: stoichiomet- 

ric H 2 -O 2 at pi = 1 atm, 7', = 300 K. p a = 0.1 atm (0 = 10). 

Figure 7: Effect of blowdown pressure ratio on specific impulse for PDRE geometries A 

and B. Propellant mixture: stoichiometric H 2 -O 2 at p, = 1 atm, 7/ = 300 K. 

Figure 8: Effect of C-D nozzle area expansion ratio on specific impulse for PDRE geometry 

C (upper panel) and D (lower panel). The dashed line in each graph represents the fixed e 
giving optimum performance at each 0. Propellant mixture: stoichiometric H 2 -O 2 at Pi = 
1 atm, Ti = 300 K. 

Figure 9: Specific impulse plotted as a function of volume ratio. Upper panel: 0 = 1000. 

Lower panel: 0 = 10. Propellant mixture: stoichiometric H 2 -O 2 at pi = 1 atm, Ti = 300 K. 

Figure 10: Normalized blowdown time (tbdUc-j / Ldt) plotted as a function of blowdown pres- 
sure ratio for PDRE geometries A and B. Propellant mixture: stoichiometric H 2 -O 2 at pi = 
1 atm, Ti = 300 K. 

Figure 11: Normalized time ( Wc-j/Ldt ) plotted as a function of blowdown pressure ratio 
for PDRE geometries C and D. Propellant mixture: stoichiometric H 2 -O 2 at Pi = 1 atm, T,- = 
300 K. 

Figure 12: Schematic of SSRE geometry studied. 

Figure 13: Performance comparison of various pulse detonation and steady-state rocket 
devices, using finite-rate quasi-l-D CFD calculations. PDRE performance is based on a single 
blowdown shot. Propellant mixture: stoichiometric H 2 -O 2 at Pi = 1 atm, 7/ = 300 K. 
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